Electronic structure of A-type antiferomagnetic LaMnOs by GW approximation 



Y. Nohara 1 , A. Yamasaki 2 , S. Kobayashi 3 , and T. Fujiwara 1,4 
1 Department of Applied Physics, University of Tokyo, Tokyo 113-8656, Japan 
2 Max-Planck Institut fur Festkorperforschung, Heisenbergstrasse 1, Stuttgart D-70569, Germany 
3 Texas Instruments Japan Ltd., 6-24-1 Nishi-Shinjuku, Shinjuku-ku, Tokyo 160-8366, Japan 
4 Core Research for Evolutional Science and Technology ( CREST-JST), Japan Science and Technology Agency 

(Dated: February 2, 2008) 

The electronic structure of the A-type antiferromagnetic insulator LaMnC>3 is investigated by 
the GW approximation. The band gap and spectrum are in a good agreement with experimental 
observation. The on-site/off-site dynamical screened Coulomb interaction on Mn sites is also calcu- 
lated. Comparison with the results by the Hartree-Fock and the static COHSEX approximations is 
discussed in detail. The dynamical screening effect is important for the band gap and widths. 
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Transition metal oxides show particular character- 
istics of physical properties, especially large control- 
lable change of properties with a small change of elec- 
tron/hole doping or the external field. Among them, 
the hole-doped manganese oxides of perovskite structure, 
Ai^B^MnOa (A=La 3 +, Pr 3+ , Nd 3+ , • • • and B=Ca 2 +, 
Sr 2+ , Ba 2+ , • • • ), have been paid much attention because 
of their properties, such as colossal magneto resistance 
(CMR), orbital, spin and charge order, metal-insulator 
transition and others. Q] 

The starting material LaMn03 is A-type antiferromag- 
netic (A-AF) at low temperatures. The Mn 3+ ion has 
t^e* configuration, whose local spin quantum number 
S=2 caused by the Hund coupling. The atomic spin cou- 
ples ferromagnetically with those of neighboring Mn ions 
on a same basal plane and couples antiferromagnetically 
with those of Mn ions on different planes. 0] The system 
shows the C-type orbital ordering with the Jahn- Teller 
(JT) distortion. Q 

The local spin density approximation (LSDA) is pow- 
erful tool to investigate the electronic structures. How- 
ever, the LSDA sometimes does not give an overall pic- 
ture of the electronic structures in transition metal com- 
pounds with strong correlation. In fact, the electronic 
structures of LaMnC>3 were calculated by the LSDA and 
the LSDA+U methods and these methods reveal sev- 
eral problems. For example, the LSDA underestimates 
the band gap and overestimates d-band width. Q In 
the LSDA+U method, [lH] the Hubbard-type Coulomb 
interaction U is introduced so that the resultant band 
gap agrees with the observed optical gap. However, this 
method produces a strong bonding state with O-p and 
Mn-d orbitals in the deep valence region as an artifact. 

The GW approximation (GWA) is a promising method 
for systems with electron-electron correlation but re- 
quires much heavier computational task. The GWA 
was applied to LaMn03 with a smaller unit cell in fic- 
titious cubic paramagnetic phase, neglecting the JT dis- 
tortion. However, the JT distortion and tilted MnOe 
octahedra (GdFeGvtype distortion) are crucially impor- 
tant in A-AF LaMnC>3. || The electronic structure of 
A-AF LaMn03 has not been calculated so far, because 



the unit cell contains 20 atoms, i.e. 4 LaMn03, and large 
memory space and the parallel computation arc needed. 

In this letter, the electronic structure of the A-AF 
LaMnC"3 is investigated by the GWA, adopting the exper- 
imentally determined crystal structure of the orthorhom- 
bic Pbnm symmetry with the JT and the GdFe03-type 
distortion. [!j Present calculation is based on the Lin- 
car Muffin-Tin Orbital method with the Atomic Sphere 
Approximation (LMTO-ASA). [ljj We use the 4 x 4 x 4 k- 
mesh in the Brillouin zone and the resultant 30 k-points 
in the irreducible zone. The set of the maximum orbital 
angular momentum of the LMTO basis in La, Mn, O and 
empty spheres are chosen to be (fddp). All E v 's for the 
LMTO linearization energy points are fixed at the av- 
eraged energies of occupied bands for respective orbital 
except La-p, for which E v is chosen in the unoccupied 
region. 

The GWA is based on the many-body perturbation 
theory and is actually the first term approximation of 
the perturbation series for the one-body Green's func- 
tion, ^3 where the self-energy is replaced by the lowest 
order term of the expansion as £ = iGW . Here G and 
W are the one-particle Green's function and the dynami- 
cally screened interaction, respectively The dynamically 
screened interaction W in the GWA is treated by the 
random phase approximation (RPA) as 

W = v + v X °W, (1) 

where v is the bare Coulomb interaction and x° is the 
irreducible polarization function x° = —iGG. 

We presume the wave- functions {V'kn(r)} of the LSDA 
to be reasonably good starting wave-functions, and then 
we adopt the LSDA Hamiltonian H = T + F H + V^ DA 
as an unperturbed one. Here T is the kinetic en- 
ergy, V H is the Hartree potential, and V^g DA is the 
exchange-correlation potential in the LSDA. The self- 
energy correction can be written with three terms as 
AE = S x + S c - V^ DA , where E x (= iGv) is the exchange 
part (the Fock term) and E c (= iGW c ) is the dynamical 
correlation part. W c is the second term in eq.QJ. The 
Green's function is defined as 



G=(E-H - AE)" 



(2) 
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and the resultant quasiparticle energy is given as 

ZknReAEk n (ekn), 



(3) 



where e^ n is the LSDA eigen-energy. The self-energy 
is AE k „(e kn ) = (V-k„|S x + S c (e k „) - ^SDA^kn) 
and the renormalization factor is Z^ n = (1 — 
9ReASi C n(e k n)/5£') _1 . The renormalization factor Z^ n 
is an inverse effective mass ratio. In the present 
work we perform one-iteration calculation without self- 
consistency, i.e. £ = iGqWq and x° = —iGoGo- 
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FIG. 1: Energy band of the A-AF LaMnOs by the LSDA 
(a), and by the GWA (b). The energy zero-th is fixed at the 
top of the respective valence band. 

The important steps of our GWA are the offset 
method, the product-basis set ^| and parallel com- 
putation. The offset method is essential to reproduce 
a response function with a small momentum transfer Kk 
and we use nearby 6 k offset points ^(±0.0784 x f , 0, 0), 
^(0, ±0.0784, 0), ^(0, 0, ±0.0784 x°f ) instead of k = 0, 
where a, b and c are the lattice constants in the or- 
thorhombic system. In the calculation of the self-energy, 
a product of two spherical wave-functions 4>4> is expanded 
in terms of the new spherical functions called a prod- 
uct basis, whose total number is 1240 for the present 
choice. The product bases are specified by the total an- 
gular momentum of two wave- functions and we adopt the 
maximum total angular momentum of the product bases 
up to (fddp) for the calculation of S c as the LMTO's 
themselves. This choice of total angular momentum re- 
duces the number of the product bases to 700. Sec also 
refs-HHH. 

The band structures calculated by (a) the LDA (ekn) 
and (b) the GWA (E kn ) are shown in Figure^ The band 
gap by the LSDA is as small as 0.16 eV and, on the other 
hand, that by the GWA becomes 0.82 eV, which agrees 
reasonably well with experimental optical gap 1.1 eV. [lhl ] 
We can see that the width of the majority Mn-d(t2 g ) 
band (in the energy region —2.3 eV < ekn < —1 eV 
in the LSDA results) becomes narrower, and that of 
the majority Mn-d3 Z 2_ r 2 band (in the energy region 
— 1.0 eV < ek„ < eV in the LSDA results) becomes 
broader. The La-d and La-f bands appear in the LSDA 
calculation at the energy range of 3 eV < ekn < 7.5 eV 
above the top of the valence band. In the GWA, these 



bands, La-d and La-f, are very strongly affected and shift 
to the energy range of 10 eV < e kn < 15 eV. For the La-f, 
the band positions in the LSDA seem too low and also the 
LSDA Hamiltonian might not be a good starting point. 
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FIG. 2: The renormalized self-energy Zk n ReA£kn(ekn) as a 
function of the LSDA eigen-energies ekn. The energy zero-th 
is fixed at the top of the LSDA valence band. 



The renormalized self-energy ZHeAT, is shown in Fig- 
ure |2l as a function of the LSDA eigen-energy ekn ■ A 
step can be seen at ekn = and it makes the band 
gap wider. In an energy region where ZRcAS is de- 
creasing/increasing as a function of £k n , the width of 
these bands becomes narrower/broader. Thus, we find 
that, in the GWA, the O-p band in the energy range 
-7.5 eV < e kn < -2.5 eV and the Mn-d(t 2g ) band in the 
energy range —2.5 eV < e kn < — 1.0 eV become narrower. 
The reduction of the O-p band is attributed to that of 
hybridized Mn-t2 g component. The value of Z is around 
0.6 in the Mn-d band region (-2.5 eV < E < 3.0 eV) 
and around 0.7 in the region of the O-p valence bands. 
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FIG. 3: (Color online) (a) Total DOS of the LSDA, (b) total 
DOS of the GW quasiparticles, and the imaginary part of 
the Green's function — (l/7r)ImG for (c) the total and (d) 
the partial components per one respective atom. In (d), the 
upper and the lower panels show the partial components with 
majority and minority spins, and bold solid line: Mn-d 3 ,!_ r 2, 
thin solid line: Mn-d 1 2_ J 2, bold dotted line: Mn-d(t2 9 ), thin 
dotted line: O-p, bold chain line: La-d, dot-chain line: La-f. 
The enrgy zero-th in (a) is fixed at the top of the valence 
bands of the LSDA and in (b)~(d) at that of the GWA. 
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Figure[3]shows the density of states (DOS) or the imag- 
inary part of the Green's function by the LSDA and the 
GWA. The band gap in the total DOS for the quasipar- 
ticle energy (b) is 0.82 eV. The imaginary part of the 
Green's function shows the smoothed and widely spread 
spectrum, because a quasiparticle has a finite lifetime. It 
should be noticed that the imaginary part of the Green's 
function gives an excellent agreement, except the La-d, 
with the observed XPS and XAS. {[ft 

In Figure^d), the notation of orbitals refers to the lo- 
cal coordinate. For example, d 3z 2_ r 2 is the orbital along 
the longest Mn-0 pair axis and d x 2_ y 2 on the perpen- 
dicular plane nearly along other Mn-0 pair axes. The 
orbitals of d(t2 S ) and d(e fl ) are split by the cubic crys- 
tal field and those of d^ z 2_ r 2 and A x 2_ y 2 mainly by the 
JT distortion. The width for the valence bands does not 
change much and, on the other hand, that of the conduc- 
tion bands becomes much broader. The partial Green's 
function is identical for all Mn ions, once one uses local 
orbitals in accordance with the local coordinates. This 
facts clearly show the characteristic orbital order (orbital 
C-type 'antiferro'-coupling) in LaMn0 3 ; i.e. each local 
orbital d 3z 2_ r 2 aligns along the elongated Mn-0 axis with 
antifcrro- and ferro-coupling within and between basal 
planes, respectively. 
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FIG. 4: The self-energy (left) and Green's function (right) 
of the highest valence (top) and lowest conduction (bottom) 
bands at k = (solid line: real part, dotted line: imaginary 
part). The energy zero-th is set at the top of the LSDA va- 
lence band. The inset shows the self-energy in a narrower en- 
ergy range and the arrows indicate the position of the LSDA 
eigen-energies. 



Figure 0] shows the self-cncrgics and the Green's func- 
tions of the top of the valence band (d 3z 2_ r 2) and the 
bottom of the conduction band (majority d x 2_ y 2 + mi- 
nority d(t2 S )) at T-point (k = 0). Actually the overall 
feature of the self-energy is not very sensitive on the k- 
points. One can notice the large difference of the shape 
of the Green's function; the one of the conduction band is 
much broader than that of the valence band. This fact is 
also confirmed for all deeper valence bands. The broad- 
ness of the Green's functions of the conduction bands is 
consistent to the behavior of the imaginary part of the 
self-energy. This fact is presumably attributed to the 
fact that the wave-functions of the conduction bands ex- 
tent appreciably over the neighboring oxygen sites and 
are affected by strong polarization fluctuation. 

The magnitude of the on-site/off-site d-d Coulomb in- 
teraction on Mn sites (<j>d<pd\W (E)\(f>d4>d) is shown in Fig- 
ure El The off-site one is for a nearest- neighbor Mn pair 
on a basal plane (ferromagnetically coupling pair) and 
the other one for a nearest-neighbor Mn pair on differ- 
ent planes (antifcrromagnctically coupling pair) is almost 
identical. The Coulomb interaction is strongly screened 
at low energies and is reduced to a few eV; 3.00 eV 
for on-site one with majority spins, 2.40 cV for on-site 
one with minority spins, 0.34 eV for off-site ones both 
with majority and minority spins. On the other hand, 
it approaches at the high energies to the values of bare 
ones 23.4 eV, 17.2 eV, 3.7 eV, respectively, because of no 
screening mechanism at high energy region. 
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FIG. 5: The dynamical screened Coulomb interaction in 
LaMnOs (solid line: interaction with majority spins, dotted 
line: interaction with manority spins). The arrows indicate 
the extremum values at the higher energies. 

The screening effects depend upon the excitation en- 
ergy, since the electrons participating in the screening 
mechanism depend very much on their excitation energy. 
The screening for the on-site Coulomb interaction is quite 
efficient and this phenomena is attributed to the screen- 
ing by easily mobile e s electrons. @ The d(t2 S ) wave- 
functions is more localized than the LSDA picture and 
the width of the Mn-d(t2 3 ) bands becomes narrower as 
seen in Figure ^b). The Coulomb interaction with mi- 
nority spins is smaller than that with majority spins, be- 
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cause the radial wave-function 4'a{ t ) °f minority spin has 
much components away from the core region and this 
fact results in a smaller values of a bare on-site Coulomb 
integral. It should be noticed that the screened Coulomb 
interaction is reduced to W ~ around E ~ 10 eV in 
Figure E| The excitation from O-p occupied states to La 
levels (3 ~ 7.5 eV) hybridized with Mn-d causes strong 
screening effects for Mn-d electrons. However, this large 
drop of the Coulomb interaction may be problematic, be- 
cause the position of La-f bands may be too low in the 
LSDA calculation. Presumably, the large reduction of W 
in the energy region of La-f or La-d bands could be true 
but it should happen at higher energy region. 

TABLE I: The band gap Ea, the band widths of occupied 
d 3z 2_ r 2, d(t2 B ) and O-p (W 3z 2_ r 2, Wt 2g and Wo p ) and the 
magnetic moment M(/j,b)- The band gap and widths are ob- 
tained from the calculated quasiparticle energy. The magnetic 
moments by the HFA, the COHSEX and the GWA are all the 
same. 
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GWA 
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E a (eV) 


0.16 


10.04 


0.71 


0.82 


1.1° 


W 3z 2_ r2 (eV) 


0.9 


2.8 
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1.1 




W t2g (eV) 
Wo P (eV) 


1.1 


1.3 


1.0 
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4.8 


5.4 


5.4 


4.6 




MGub) 


3.49 








3.7\ 3.87 c 



°Ref. E3. 

'j.B.A.A. Elemans, B. van Laar, K.R. van der Veen and B.O. 
Loopstraj J. Solid State Chem. 3, 238 (1971). 
c Ref. 0. 



We summarize in Table 1 the calculated results for 
band gap, widths of occupied bands and magnetic mo- 
ment with experimental results for comparison. The cal- 
culated results are for the LSDA, the Hartree-Fock ap- 
proximation (HFA), the static COHSEX approximation 
(sCSA) [13 and the GWA. The sCSA is obtained with 
the self-energy of E = or the zero excitation energy . 
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The band gap is found to be 0.82 eV by the GWA, which 
is in good agreement with experimental value 1.1 eV and 
slightly larger than the value by sCSA. The large band 
gap of the HFA is owing to the approximation free from 
the self-interaction. Then the sCSA adds the effect of 
the static screening. In fact, the change from the HFA 
to the cCSA is very large for the unoccupied Mn-d band. 
The dynamical part of the electron-electron interaction 
is attributed to the small difference between the band 
gaps of the sCSA and the GWA. The magnetic moment 
of magnetic insulators is not changed from the results of 
the LSDA in the HFA, the sCSA and the GWA, because 
these methods are the perturbative ones and we do not 
renew wave-functions. The band widths of Mn-d 3z 2_ r 2 
states by the GWA and the sCSA are slightly greater 
than that of the LSDA. The band width of Mn-d(t 2ff ) 
states by the GWA and the sCSA are less than that of 
the LSDA. The band width of O-p bands by GWA is 
less than the LSDA value. These results of the band 
widths indicate the wave- functions of Mn-d(t2 ff ) and O-p 
states becomes more localized by the GWA and that of 
Mn-d(e g ) less. Therefore, the screening effect affects to 
d(t2 S ) and d(e fl ) states in an opposite way. 

In summary, we applied the GWA to the A-type anti- 
ferromagnetic LaMnOa. Band gap and magnetic moment 
are in good agreement with experimental results. The 
band width of Mn majority d(t 2s ) states become nar- 
rower due to the correlation effect and, on the contrary, 
that of Mn majority d3 Z 2_ r 2 states becomes broader. 
The spectrum shows an excellent agreement with exper- 
imental results. We also present the energy depending 
Coulomb interaction, and hnd that on-site d-d Coulomb 
interaction is strongly screened at the low energy region 
due to mobile e g electrons. 

Computation was partially carried out by use of the 
facilities at the Supercomputer Center, Institute for Solid 
State Physics, University of Tokyo, and the Institute of 
Molecular Science at Okazaki. 
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